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Abstract. Neutrino-pair bremsstrahlung by relativistic 
degenerate electrons in a neutron-star crust at densities 

g cm~ 3 is analyzed. The 



10 9 g cm ~ 3 < p & 1.5 x 10 1 
processes taken into account are neutrino emission due 
to Coulomb scattering of electrons by atomic nuclei in 
a Coulomb liquid, and electron-phonon scattering (the 
phonon contribution) and Bragg diffraction (the static- 
lattice contribution) in a Coulomb crystal. The static- 
lattice contribution is calculated including the electron 
band-structure effects for cubic Coulomb crystals of dif- 
ferent types and also for the liquid crystal phases com- 
posed of rod- and plate-like nuclei near the bottom of 
the neutron-star crust (10 14 g cm~ 3 ^ p ^ 1.5 x 10 14 g 
cm -3 ). The phonon contribution is evaluated with proper 
treatment of the multi-phonon processes which removes 
a jump in the neutrino bremsstrahlung emissivity at 
the melting point obtained in previous works. Generally, 
bremsstrahlung in the solid phase does not differ signif- 
icantly from that in the liquid. At p <, 10 13 g cm~ 3 , 
the results are rather insensitive to the nuclear form fac- 
tor, but results for the solid state near the melting point 
are affected significantly by the Debye- Waller factor and 
multi-phonon processes. At higher densities the Debye- 
Waller factor and multi-phonon processes become less im- 
portant but the nuclear form factor becomes more sig- 
nificant. With growing p, the phonon contribution be- 
comes smaller. Near the bottom of the neutron star crust 
bremsstrahlung becomes less efficient due to the reduc- 
tion of the effective electron-nucleus matrix element by 
the electron band-structure effects and the nuclear form 
factor. A comparison of the various neutrino generation 
mechanisms in neutron star crusts shows that electron 
bremsstrahlung is among the most important ones. 
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1. Introduction 

Neutrino-pair bremsstrahlung of electrons in liquid and 
crystalline phases of dense matter is one of the major neu- 
trino energy-loss mechanisms in neutron star crusts. Here, 
by bremsstrahlung we mean neutrino emission due to elec- 
tromagnetic interaction of electrons with atomic nuclei. 
The process can be written schematically as 



(Z,A) 



(Z,A) 



(1) 



ft proceeds via neutral and charged electroweak currents 
and leads to emission of neutrinos of all flavors. 

For practical application to the thermal evolution of 
neutron stars, one needs to know the neutrino energy emis- 
sion rate (emissivity) Q in the density range from about 
10 9 g cm~ 3 to 1.5 x 10 14 g cm -3 (the core-crust interface) 
at temperatures T 5 x 10 9 K (at which the nuclei are 
not dissociated). Under these conditions, the electrons are 
strongly degenerate and ultra-relativistic, and the nuclei 
form either a Coulomb liquid, or a Coulomb crystal. For 
densities higher than 10 12 - 10 13 g cm -3 , the melting tem- 
perature of the crystal exceeds 5 x 10 9 K, and the case of a 
Coulomb liquid is of no practical importance. In the den- 
sity range from about 10 14 g cm -3 to 1.5 x 10 14 g cm" 3 , 
the nuclei resemble rods and plates, rather than spheres 
(Lorenz et al. |1993| ; Pethick & Ravenhall |l995| ). 

The neutrino-pair bremsstrahlung process (Q) in a 
crystal is formally different from that in a liquid. In the 
liquid state, neutrinos are generated due to Coulomb scat- 
tering of electrons by nuclei. In the solid state, there are 
two contributions to the process, electron-phonon scatter- 
ing (electron scattering by the nuclear charge fluctuations 
due to lattice vibrations, referred to as the phonon con- 
tribution), and the Bragg diffraction of electrons, which is 
commonly called the static-lattice contribution. 

Neutrino-pair bremsstrahlung has been analy z ed by 
a number of authors (see, e.g., Itoh et al. 1989 , 1996 ; 
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Pethick & Thorsson |1997| , and references therein). The 
case of a Coulomb liquid has been thoroughly studied 
by Festa & Ruderman ( J1969Q , Dicus etal. ( p. 9761 ), Soyeur 
& Brown ( p.979| ), Itoh & Kohy ama ( |1983| ), and, most re- 
cently, by Haensel et al. ( |1996| ). The phonon contribution 
in the crystalline lattice has been analyzed by Flowers 
( |1973D , Itoh e t al. ( |l984b| , |1989[ ), and also by Yakovlev 
& Kaminker (199£). Notice, however, that all these au- 



thors have made use of the so-called one-phonon approx- 
imation, whereas, as we shall demonstrate below, multi- 
phonon scattering processes can give a noticeable contri- 
bution especially near the melting point. The static -lattice 
contributio n has b een considered by Flowers ( 1973 ) and by 
Itoh et al. (|l984a ) neglecting the electron band-structure 
effects (the presence of gaps in the elect ron d i spersi on re- 
lation). However, Pethick & Thorsson ( |l994 |1997| ) have 
shown that the band structure can suppress strongly the 
static-lattice contribution and is thus very important. In 
particular, Pethick & Thorsson (1997) derived a general 
expression for the static-lattice contribution. 

In the present article, we give an overall analysis of 
neutrino-pair bremsstrahlung by electrons in neutron star 
crusts, and obtain formulae convenient for applications. 
In Sect. 2 we describe briefly physical conditions in neu- 
tron star crusts. In Sect. 3 we summarize the formalism 
for calculating neutrino-pair bremsstrahlung in the liq- 
uid and solid phases. While studying the solid phase we 
consider the ordinary body-centered-cubic (bec) crystals, 
as well as face-centered-cubic (fee) and hexagonal-close- 
packed (hep) cubic crystals which can also be formed in 
neutron-star crusts. The phonon contribution will be cal- 
culated using a new approach (Baiko et al. 



1998) 



incor- 
porating multi-phonon effects. It eliminates a jump in the 
neutrino emissivity at the melting point obtained in pre- 
vious articles. The static-lattice contribution will be an- 
alyzed not only for the traditional phase of spherical nu- 
clei, but also for the 'exotic' phases of nonspherical nuclei. 
In Sec. 4 we examine the main properties of neutrino- 
pair bremsstrahlung at various densities and temperatures 
in neutron-star crusts, and compare bremsstrahlung with 
other neutrino emission mechanisms. We present also a 
simple analytic fit for practical evaluation of the neutrino 
bremsstrahlung energy loss rate. 



2. Physical conditions 

Let us outline physical conditions in a neutron star enve- 
lope at densities 10 9 g cm -3 ^ p ^ 1.5 X 10 14 g cm" 3 , 
and at temperatures T ^ 5 x 10 9 K (Sect. 1). Matter 
in these layers consists of electrons and atomic nuclei 
(ions). At densities higher than the neutron drip density, 
Pd ~ 4.3 x 10 11 g cm -3 , free neutrons appear between 



the nuclei (e.g., Negele & Vautherin |1973j ). At p ^ 10 14 g 
cm" , the nuclei are likely to form nonspherical clusters 



(Lore nz et al. 1993 Oyamatsu 1993 ; Pethick & Ravenhall 
1991) . 



The state of degenerate electrons is characterized by 
the Fermi momentum p-p = Tik-p or the relativistic param- 
eter x: 



Pl- 



pp = ?i(37rV) /J , x = — « 100.9 (p 12 Y c ) 



1/3 



(2) 



where Y c — ric/n^ is the number of electrons per baryon, 
n e is the number density of electrons, rib the number den- 
sity of baryons, and p\2 is mass density in units of 10 12 g 
cm" 3 . In the density range under study the electrons are 
ultra- relativistic (x ^S> 1). The electron degeneracy tem- 
perature is 



TV 



(Vl+x 2 -l)T , T = 



m c c 



5.930 xlO 9 K, (3) 



where fee is the Boltzmann constant. In our case T ^ 
5 x 10 9 K, and the electrons are strongly degenerate. 

The nuclear composition of neutron-star envelopes is 
not very well known, although it is quite certain that light 
elements such as H and He transform into heavier ones at 
densities which are lower than the densities of interest. 
For simplicity, we assume that only one nuclear species 
is present at any fixed density (pressure). This leads to 
discontinuous variations of the nuclear composition with 
density (pressure). The temperature dependence of the nu- 
clear composi tion c an be ignored at T < 5 x 10 9 K (e.g., 
Haensel et al. 1996). For illustration, we shall make use of 



two models of matter in a neutron-star crust: ground-state 
(cold-catalyzed) matter and accreted matter. For describ- 
ing ground-state matter, we sha ll use the following data: 
the results of Haensel & Pichon ( 1994 ) at p < pd based on 
new laboratory measurements of nuclear masses wit h large 
neutron excess, the results of Negele & Vautherin ( 1973 ) 
for spherical nuclei at pd ^ p ^ 1.5 x 10 14 g cm" 3 de- 
rived by a modified Hartree-Fock method, and model I of 
Oyamatsu ( |1993| ), which takes into account nonspherical 
nuclei. We shall adopt the composit ion of accreted mat- 
ter calculated by Haensel & Zdunik (199C) for p <, 10 13 g 



cm , and we shall not consider accreted matter at higher 
densities. Accreted matter consists of lighter nuclei, and 
neutron drip is shifted to higher density as compared to 
ground-state matter (Fig. |l|). 

The state of spherical nuclei is determined by the ion- 
coupling parameter 



Z 2 e 2 Z 5 / 3 
T = — — w 0.225a;- 



r 8 



(4) 



where Ze is the nuclear charge, a = [3/(47rni)] 1 / 3 is the 
ion-sphere (Wigner-Seitz cell) radius, rt; = n c /Z is the 
number density of nuclei, and Tg is the temperature in 
units of 10 8 K. For the densities and temperatures of in- 
terest, the spherical nuclei constitute either a strongly- 
coupled Coulomb liquid (1 < T < T m ), or a Coulomb 
crystal (r > r m ), where r m =172 corresponds to solidifi- 
cation of a classical one-component Coulomb liquid into 
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a bcc lattice (Nagara et al. 1987). Thus the melting tem- 
perature is 



T m = 



Z 2 e 2 
afcer E 



1.32 x 10 7 Z 5/3 (pi 2 Y c ) 1/3 K 



(5) 



The density profiles of T m for the ground-state and ac- 
creted matter are presented in Fig. [j]. The melting tem- 
perature of accreted matter is systematically lower due 
to the lower values of Z. If p > 10 13 g cm" 3 , one has 
T m <; 5 x 10 9 K for the ground-state model, and matter 
is always solid for the conditions under study. A classi- 
cal bcc lattice is bound most tightly (see, e.g., Brush et 
al. 1966). Therefore, it is widely assumed that neutron- 



star crusts are composed of such crystals. However, an fee 
or hep crystal is bound only slightly more weakly, so fee 
and hep crystals may well occur in dense stellar matter 



along with the bcc ones (e.g., DeWitt et al. 1993; Baiko & 



Yakovlev 1995), and we shall also consider this possibility. 
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(g) = ground-state matter 
(a) = accreted matter 



fcuO 




lg p [g cm" ] 

Fig. 1. Density dependence of the melting temperature 
T m and the ion plasma temperature T p , and also neutron 
drip density (vertical dotted lines) for ground-state matter 
and accreted matter 

An important parameter for a Coulomb crystal is the 
ion plasma temperature 



rp P 

P_ k B 



7.83 x 10 y 



ZY c p 12 



1/2 



K. 



(6) 



which is determined by the ion plasma frequency lu p = 

1/2 

(47r.Z 2 e 2 7ii/mi) , where ra- x f=a A\m n is the ion mass and 
m u = 1.66 x 10~ 24 g is the atomic mass unit. The plasma 
temperature characterizes ion vibrations. If T <; T p /8 the 



vibrations can be treated classically (the high-temperature 
case) while at T <C T p they should be treated quantum- 
mcchanically (the low-temperature case). In neutron-star 
crusts, a transition from the high-temperature to the low- 
temperature regime takes place at temperatures which are 
lower than T m (Fig. 0). 

We shall consider ion vibrations in a purely Coulomb 
lattice. For p > Pd there are neutrons in the space between 
ions (nuclei) . How to include the effect of these neutrons 
on the lattice dynamics is at present an unsolved problem, 
and we shall ignore this effect. 

The effective interaction between electrons and nuclei 
depends on a number of effects. One is the character of 
the Coulomb interaction itself, a second is screening of 
the interaction by electrons, a third is the shape of the 
nuclear charge distribution, and a fourth is the effect of 
thermal vibrations effectively smearing out the nuclear 
charge, an effect taken into account by inclusion of the 
Debye-Waller factor. Thus the Fourier transform of the 
effective electron-ion interaction may be written as 



_ Anepz F(q) w[q) 



(7) 



where pz is the ion charge per unit volume, F(q) is the 
form factor which reflects the charge distribution within 
one nucleus, e(q) is the static longitudinal dielectric factor 



(Jancovici 1962] ), and W(q) is the Debye-Waller factor. 
For a crystal of spherical nuclei, we have pz = Zen\ ; Eq. 
(0) is valid also for non-spherical nuclei. The Debye-Waller 
factor in a crystal of spherical nuclei can be written as 



x*rt \ h 1 2 /coth(^ s /2fc B T) 
W ® = 4^\ £ 



(8) 



where uj s is a phonon frequency, s enumerates phonon 
modes and brackets (. . .) denote averaging over phonon 
wave vectors and polarizations (e.g., Eq. (8) in Yakovlev 
& Kaminker 1996 ). In a Coulomb crys tal W = W(q) is 
accurately fitted by Baiko & Yakovlev ( [1995 ) : 



W 



2 \2k F J 



o 



a 



u-xe 9 ' ltp +t p u-2 



where t p = T/T p and 



tv 



4m 2 c 2 
fc B Tpr7ii 



1.683 



AiZ 



(9) 



(10) 



The quantities U—i and U-i are the frequency moments 
of the phonon spectrum, u n = ((uj s /uj p ) n }. For the bcc 
lattice, the frequency moments are well known (see, e.g., 
Pollock & Hansen 1973). For the fee and hep crystals, 



they are easily derived by calculating the phonon spec- 
trum from the Ewald transformation of the lattice sums. 
The phonon spectra of the bcc, fee, and hep lattices ap- 
pear to be quite similar. Accordingly the properties of 
these crystals are very similar. In particular, U—% = 2.798, 
2.720 and 2.703; u_ 2 = 12.972, 12.143 and 12.015, for 
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bcc, fee and hep lattices, respectively. Note that the fre- 
quency moments for the fee lattice calculated by Baiko & 



Yakovlev (1995) and used by Baiko & Yakovlev (199£) and 
Yakovlev & Kaminker (1996) are inaccurate due to erro- 



neous boundaries of the Brillouin zone for this lattice im- 
plemented in the momentum averaging scheme. Improved 
calculations show that Eq. (|9|) remains valid with the ac- 
curate frequency moments for all the lattice types under 
study. The Debye- Waller factor is important if ion vibra- 
tions (either thermal or zero-point ones) are strong. 

The simplest nuclear form factor in Eq. (Q) corresponds 
to a spherical atomic nucleus with a uniform proton core 
of radius R c : 



F(q) 



(qRc) 3 



[sm(qR c ) - qR c cos(gi? c )] 



(11) 



At p < pd, the radii of the nuclei do not change under 
the ambient pressure of neighbouring particles, and one 
can use the standard formula R c — 1.15 A 1 / 3 fm. At these 
densities the ratio r\ = R c /a <C 1, and the effect of the 
atomic form factor is generally negligible (F(q) « 1 for 
q ^ 2fcp). At densities Pd < P ^> 10 13 g cm~ 3 the sim- 
plest form factor remains adequate, but the proton core 
radius becomes R c w 1.83 Z 1 / 3 fm as deduced by Itoh & 
Kohyama ( 1983Q from the results of Negele & Vautherin 
(1973). In this case, r\ can be as high as 0.2 - 0.3 and the 
effect of the form factor is important. For higher p, the pro- 
ton charge distribution becomes smoother, and the form 
factor should be modified. 

In order to describe accurately all layers of the neutron 
star crusts we make use of the results by Oyamatsu (1993) 
who calculated the local neutron (n) and proton (p) num- 
ber density distributions within a Wigner-Seitz cell and 
fitted them in the form 



(n 
n 



3 

out 







1 - 


(*)"] 





r<R ]7 ^ 
r>Rj, 



where j —n or p, and n™, n° , tj and Rj are the fit param- 
eters. These parameters, as well as the sizes of Wigner- 



Seitz cells, are presented in Table 6 in Oyamatsu (1993) 
for several selected values of the baryon number density 
rib f° r spherical and nonspherical nuclei. 



In particular, Oyamatsu (1993) gives the fit parame- 
ters for spherical nuclei at three values of baryon number 
density n b =0.01, 0.03 and 0.055 fm" 3 (i.e., p = 1.66 x 10 13 , 



4.98 x 10 13 and 9.13 x 10 13 g cm 3 ) in the inner neutron- 
star crust. These parameters are quite consistent with 



those presented by Negele & Vautherin (1973) for nearly 
the same rib- Some of these parameters can also be de- 
duced from Figs. 3 and 4 and from Table 3 of Negele 



& Vautherin (1973) for several other values of rib hi the 
inner crust. The parameters appear to be smooth func- 
tions of rib, so we interpolated between the given points 
at < p < 1.4 x 10 14 g cm~ 3 . This interpolation 



smears out jumps in the nuclear composition with in- 
creasing p, but these have little effect on the neutrino-pair 
bremsstrahlung. The interpolation allows us to calculate 
easily the parameters of spherical nuclei at any density in 



the inner crust. In particular, 



gives the number den- 



sity of free neutrons outside the nuclei (r > i? n ) while R n 
may be called the nuclear radius. There are no protons 



outside nuclei, n° 



0, in this regime, and the proton 



core radius R p is somewhat smaller than R n . The nuclear 
mass m; is assumed to be that of all nucleons within R n . 
The parameters t n and t p range from about 4 to 6 and de- 
crease with increasing p. Note that the proton core radius 
R p in Eq. ([l2|) is somewhat larger than the proton core 
R c in the simplified model (pi]). 

We have obtained also an analytical description of 
atomic nuclei for lower densities, 10 8 g cm -3 < p < pd 
in the outer neutron-star crust making use of the results 
by Haensel & Pichon (1994). We have adopted the same 
parameterization Jl2] ) and constructed simple analytic ex- 
pressions for the nuclear parameters of the ground-state 
matter as a function of rib in the outer crust. In this case 
n° ut = 0. At low density in the outer crust these expres- 
sions yield 56 Fe-nuclei. 

According to model I of Oyamatsu ( |l993| ), the phase 
with spherical nuclei in the inner crust is realized up to 
a density n b = 0.0586 fm" 3 (p = 0.973 x 10 14 g cm" 3 ). 
This is followed by the phase with rod-like nuclei up to 
a density n h = 0.0749 fm" 3 (p = 1.24 x 10 14 g cm" 3 ), 
and the phase with slab-like nuclei (up to rib = 0.0827 
fm -3 , p — 1.37 x 10 14 g cm -3 ). Subsequently there are 
two phases with the roles of nuclear matter and neutron 
matter reversed, the rod-like one (up to rib = 0.0854 fm~ 3 , 
p = 1.42 x 10 14 g cm -3 ), and the "Swiss cheese" (inverted- 
spheres) one, which is the analog of the phase with spher- 
ical nuclei and is the last phase in the neutron-star crust 
(up to n h = 0.0861 fm -3 , p = 1.43 x 10 14 g cm" 3 ). At 
higher density the nuclei dissolve to give the uniform mat- 
ter of the neutron star core. 



In each crystalline phase of matter the Wigner-Seitz 
cell has its own geometry, but we shall assume that in the 
case of bcc phases it may be approximated by a sphere, 
and in rod-like phases by a right circular cylinder. We 
shall assume that the nucleon density distributions may 
be described by Eq. ([l2|), where r is the distance from 
the cell center (e.g., from the axis of a rod in the rod- like 
phase, or from the symmetry plane in the slab- like case). 
We interpolate these parameters as functions of rib within 
each phase separately. In the phases with spheres, rods 
and slabs, n° ut = 0, and n° ut describes the number den- 
sity of free neutrons, and the region r < R n is occupied 
by the nucleus itself (with njf > ri° ut ). In the two last 
"bubble" phases with the roles of nuclear matter and neu- 
tron matter reversed, n° ut ^ 0, and n° ut > n™, i.e., the 
local number density of neutrons and protons increases 
with distance r from the center of the Wigner-Seitz cell. 
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With increasing density the nucleon density profiles be- 
come smoother, approaching that for uniform matter. 

Thus we have obtained a simple analytic descrip- 
tion of the neutron and proton local density profiles for 
the ground-state matter throughout the outer and in- 
ner neutron-star crusts including non-spherical phases of 
atomic nuclei. This description will be used below and 
will be referred to as the smooth composition (SC) model 
of ground-state matter. Using this model one may easily 
determine the nuclear form factor F(q) numerically by cal- 
culating the Fourier transform of n p (r). Unfortunately, vi- 
brational properties are known only for crystals composed 
of spherical nuclei, apart from recent work on the elastic 



constants (Pethick & Potekhin 1998). Thus in phases with 
non-spherical nuclei we know neither the Debye-Waller 
factor nor the phonon spectrum. 

3. General formalism 

The general expression for the neutrino emissivity Q due 
to the neutrino-pair bremsstrahlung ([!]) of relativistic de- 
generate electrons in a pl asma of spherical nuclei can be 
written as (Haensel et al. 1996 ) 



S7rG 2 Z 2 e 4 C 2 
567ft 9 c 8 



(k B T) 6 niL 



3.229 x 10 11 p 12 ZY C T 8 6 L erg s _1 cm" 



(13) 



where G F = 1.436 x 10 49 erg cm 3 is the Fermi weak cou- 
pling constant, and L is a dimensionless function to be 
determined. Furthermore, C\ = C v +C 2 A + 2(C^ + C'%), 
where CV and Ca are the vector and axial-vector con- 
stants of weak interaction, respectively. We have Cy = 
2 sin 2 9w + 0.5 and Ca = 0.5, for the emission of electron 
neutrinos, to which both charged and neutral currents con- 
tribute, and Cy = 2 sin 2 8w — 0.5 and C' A — —0.5, for the 
emission of muonic and tauonic neutrinos, which is due 
only to neutral currents. In this case, 6*w is the Weinberg 
angle, sin 2 W ~ 0.23. Equation (£|) with C\ w 1.675 
is obtained taking into account the emission of the three 
neutrino flavors {v c , v^, and v T ). 

Thus the problem reduces to evaluating the function 
L. In the liquid of atomic nuclei, L = L\ iq is determined by 
the Coulomb scattering of electrons by nuclei (Sect. 3.1). 
In the Coulomb solid, L consists of two parts, 



L — L so \ — Lph + L 8 i, 



(14) 



where L p h describes the phonon contribution (Sect. 3.2) 
and L B i describes the static-lattice contribution (Sect. 3.3). 

The non-spherical nuclei at the bottom of the neutron- 
star crust form a liquid crystal at the temperatures of in- 
terest. However, since the vibrational properties of these 
phases are largely unknown, we shall not consider the 
phonon contribution for this case. The static-lattice con- 
tribution for such a solid will be analyzed in Sect. 3.3. 



3.1. Liquid phase 

The factor Lu q is a slowly varying function of the plasma 
parameters and it has basically the same significance as 
the Coulomb logarithm in calculations of transport pro- 
cesses in plasmas. The most general expression for it was 
obtained by Haensel et al. (1996|), who found 



2A' F 



he 

xR T (q t ,q r ) R NB (q t ) 



dq r 



<? 4 K<7)| 2 



(15) 



where q = qt + q r is the momentum transfer from an elec- 
tron to a nucleus in a collision event, q t corresponding to 
purely elastic Coulomb scattering while q r takes into ac- 
count inelasticity due to the neutrino emission; F(q) is the 
nuclear form factor, e(q) is the static longitudinal dielec- 
tric function (Sec t. 2), S (q) is the ion -ion s tructure factor 
(e.g., Itoh et al. 1983| ; Young et al. 1991 ), and R^ B (qt) 
includes non-Born corrections. The function Rx{qt,qr) is 
given by Eq. (20) in Haensel et al. ( 1996 ) and describes 
the effects concerned with the thermal smearing of the 
electron distribution function. 

In the ultra- relativistic limit, the Coulomb loga- 
rithm depends actually on three dimensionless parame- 
ters: L hq = L liq (Z,-n,t F ), T] = R c /a, t F = k B T/(2p F c) « 
T/(2T F ), where T F is g iven by Eq. (|). 

Haensel et al. ( [1996 ) calculated the Coulomb logarithm 
@ with the form factor @ at Z < 50, t F < 0.1 and 
•q % 0.2 and fitted the results by an analytic formula (their 
Eq. (25)). In this article, we shall calculate Li; q from the 
starting equation ( |l5| ) since we shall not restrict ourselves 
to the simple form factor ( pi] ) (Sect. 2). We shall use the 
structure factor S{q) obtained by Rogers & DeWitt (un- 



published) and accurately fitted by Young et al. ([1991]). 

Notice that if we neglect the thermal smearing of the 
electron distribution, then Eq. ( p5|) r educe s to the familiar 
expression (e.g., Festa & Ruderman 1969) 



L\ia — 



S(q)\F(q)f 

yHq)\ 2 



2y 2 



In j/ #nb(<7), (16) 



where y = q/(2k F ). As shown by Haensel et al. (1996), it 



is a good approximation for T -C hcq s , where q s ~ a 
the Coulomb screening momentum. 



is 



The factor i?NB(<?) represents (Haensel et al. 1996) the 
ratio of the electron scattering cross sections by atomic nu- 
cleus calculated exactly and in the Born approximation. It 
describes the non-Born correction to the Born approxima- 
tion. To simplify consideration of the non-Born corrections 
we have introduced the mean non-Born correction factor 
-Rnb defined as 



r> _ rNB / r Born 
-KNB - ^li q / % q , 



(17) 



where L^® and L®° Tn are given by Eq. ( |l5] ) calculated with 
an accurate factor i?NB and with i?NB = 1 7 respectively. 
We have evaluated -Rnb using the form factor ( |TT| ) for wide 



6 



A.D. Kaminker et al.: Neutrino-pair bremsstrahlung 



ranges of the parameters 77, tp (or T) typical for neutron 
star envelopes, .Rnb appears to be a very slow function of 
77 and r. Since our treatment of the non-Born corrections 
is rather phcnomenological anyway, we have set r = 150 
and r] = 0.1 and neglected the thermal smearing of the 
electron distribution function. Then -Rnb is a function of 
the only remaining parameter, Z. The numerical results 
for Z < 60 are accurately fitted by 



(18) 



i?NB = 1 + 0.00554 Z + 0.0000737 Z 2 



This formula enables us to calculate Ln q from Eq. ( |15[ ) or 
( |l6| ) with i?NB(<?) = 1, i-e., in the Born approximation, 
and introduce the mean non-Born correction (|l8|) using 

An accurate calculation of the non-Born corrections 
in crystalline matter (Sees. 3.2 and 3.3) is a difficult task 
which goes beyond the scope of the present paper. How- 
ever, we shall see (Sect. 4) that neutrino bremsstrahlung 
in crystalline matter is quite similar to that in a Coulomb 
liquid. Thus we adopt the same factor ( |l8l ) to account for 
the non-Born corrections in Coulomb crystals. 



3.2. Phc 



contribution 



The phonon contribution in a Coulomb crystal of spheri- 
cal nuclei was studied by a number of authors (e.g., Flow- 
ers |l973j ; Itoh et al. |l984b| ; Yakovlev & Kaminker |1996| ). 
So far, all articles have been restricted to consideration 
of one-phonon processes (absorption or emission of one 
phonon). To allow for the background lattice vibrations 
the one-phonon reaction rate has usually been multiplied 
by e~ 2W , where W = W(q) is the Debye- Waller factor 
introduced in Eq. (||). 

Under astrophysical conditions at not too low temper- 
atures, the main contribution to elect ron-phonon scat- 
tering comes from umklapp processes, in which the 
electron momentum transfer Tiq in a scattering event 
lies outside the first Brillouin zone. Then the phonon 
(quasi)momentum is determined by reduction of q to the 
first Brillouin zone. The umklapp processes require q ^ qo, 
contrary to the normal processes in which q remains in the 
first Brillouin zone and q ^ qo, where qo w (67r 2 n i ) 1 / 3 is 
the radius of the Brillouin zone in the approximation in 
which it is treated as a sphere. Umklapp processes domi- 
nate since the parameter y = q /(2kF) = (4Z) -1 / 3 is typ- 
ically small (due to the large Z). Accordingly the phase 
space associated with umklapp processes is much larger 
than that for the normal ones (e.g., Raikh & Yakovlev 
1982] ) and, in most cases, it is sufficient to consider umk- 



lapp processes alone. 

Let us rederive the expression for L p h with a proper 
treatment of multi-phonon processes. We start from the 
general integral expression (Eq. (18) of Flowers 1973[ ) for 
the neutrino emissivity due to electron-phonon scattering. 
The integrand contains the inelastic part Sd(q, O) of the 
dynamical structure factor of ions in a Coulomb crystal. 



It is the inelastic part that is responsible for the electron- 
phonon scattering. Baiko et al. ( 1998j ) have obtained its 
expression for q <; go by accurate summation of multi- 
phonon diagrams: 

r+00 

S d (q,Q) = / dte mt S(q,t), 



S(q,t) = 
where 



-2W 



(e*<*> - l) 



cos [ui s (t + ih/2k B T)] 



(19) 



(20) 



uj s sinh (hu) s 1 '2fceT) 

In these calculations the density operators are calculated 
to all orders in the phonon creation and annihilation oper- 
ators, but the phonon dynamics is treated in the harmonic 
approximation. 



We now use the approach of Flowers (1973) with the 
dynamical structure factor ( |l9| ) and make the same sim- 
plifications as in deriving the expression for L p h by a 
semianalytical method described by Yakovlev & Kaminker 
( |1996| ). Then we obtain 

1 dy S eS (q)\F(q)\ 2 



'yo 

In this case 
Ses(q) 



1 



2y 2 



lny 



(21) 



63?i b 



16tt 7 (fc B T) 6 
Q + oj 



+00 



-Hoc- 



dt 



AQ,t 



S(q,t), 



(22) 



e h(n+ui)/k B T _ 1 
where hu is the neutrino-pair energy. As in Yakovlev & 



Kaminker (1996) the lower integration limit ?/o excludes 
the low-momentum transfers in which the umklapp pro- 
cesses are forbidden. 

Comparing Eqs. ( ^l| ) and ( |l6|) we see that S e s (q) plays 
the role of an effective static structure factor that defines 
the phonon contribution to the neutrino bremsstrahlung. 
The expression (^2|) for it can be easily simplified. First we 
can integrate over Q and u> which leaves us with a single 
integration over t: 

4 

„-2W 



189 



xlm 



dt 



„*(*) 



1 



(23) 

t 5 sinh 2 (7rffc B T/?i) 
The latter integration is non-trivial since the integrand 
is singular at t — 0. However, the singularity is easily 
removed by using the theory of functions of a complex 
variable. Since the function <!>(£) is analytic, the integrand 
allows us to shift the integration path into the complex 
t plane. The appropriate shift is t = t' — ih/(2h&T). It 
transforms Eq. ( p3| ) into a rapidly converging integral 

,S-.,rUy) = 189 (-) / ,l<- 



189 Ou 



1 - 40£ 2 + 80£ 4 



(l + 4£ 2 ) 5 cosh"«) 



1 



(24) 
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where £ = t'k^T/% and 

Tiq 2 I cos(w s £') 



2 ni] \ u.v sin! i ! 2k i-T) / 



(25) 



The phonon averaging (. . .) can be performed using 
the method of Mochkovitch & Hansen ( 1979 ). Afterwards 
the integral ( pif ) can be calculated numerically. The be- 
haviour of S c g(q) depends on temperature, and may be 
characterized by the dimensionless parameter t p — T/T p 
introduced in Eq. (0). First consider the asymptotes for 
t p <C 1 and t p <; 1. 

In the low-temperature case, t p <C 1, it is sufficient 
to set e* - 1 = $ in Eqs. ©, (||) or (||). This case 
corresponds to the familiar one-phonon approximation 
adopted in previous articles. In the notation introduced 
by Yakovlev & Kaminker (1996) one finds in this case 
that 



Sl§ (q) = a y 2 t p G(t p ) e" 



2W 



a y 2 bt 2 e~ 



2W 



(26) 



where ao is defined in Eq. (|10|), G(t p ) is the function deter- 
mined by Eq. (12) in Yakovlev & Kaminker (1996), and 
b is a numerical factor specified by the low-temperature 
asymptote G(t p ) = bt p . For a bcc Coulomb crystal, one 
has b w 231. This value is more accurate than b w 202 
given by the numerical fit (15) in Yakovlev & Kaminker 
(1996) since the latter authors did not intend to pro- 
duce a fitting formula which would be highly accurate for 
t v ^ 0.01. 

Therefore, the results obtained in the previous articles 
are strictly valid for T <C T p . The Debye- Waller exponent 
e -2W j ncm ded into the one-phonon approximation takes 
into account renormalization of the one-phonon interac- 
tion by background lattice vibrations. Notice that recent 
calculations of G(t p ) in Yakovlev & Kaminker ( 1996 ) for 
the fee lattice are inaccurate due to an error in specifying 
the boundaries of the first Brillouin zone (Sect. 2). After 
correcting this error, we obtain a result for G(t p ) very close 
to that for the bcc lattice. Note also a misprint in Eq. (18) 
for the function F\ in Yakovlev & Kaminker ( |1996 ): in the 
second term (1.5 + a) 3 / 4 is erroneously printed instead of 
(1 + a) 3 / 4 although all calculations were made using the 
correct expression. 

In the opposite case of high temperatures, t p PS 1, the 
asymptotic form of Eq. (124) is very simple, 



Scff(g) = 1 - e 



-2W 



(27) 



Thus S e g(q) becomes noticeably larger than in the one- 
phonon approximation, as a result of multi-phonon pro- 
cesses. 

We have also calculated S c s(q) for a wide range of tem- 
peratures T and density parameters ao defined in Eq. (|l(J) . 
We have verified that, under the conditions in a neutron- 
star crust, «o < 0.2. For such ao and all t p the numerical 
results are fitted by a simple expression 



where 
W x = 



a y 2 bu^ 2 tl 



2 J (btp) 



- 2 exp(-7.6t p ) 



(29) 



and b « 231, for a bcc lattice. Notice that W\ « W ~ 
0.5 aoy 2 u-2t p for t p ^ 1. Notice also that S e g(q) given by 
Eq. (28) reproduces quite accurately Slg h (q) if we replace 
e 2Wl — 1 by 2W%. The effect of multi-phonon processes 
on neutrino bremsstrahlung emission will be described in 
Sect. 4. 




y 



Fig. 2. The effective structure factor S e s{q) which en- 
ters the normalized neutrino emissivity ( pl| ) , produced by 
the electron-phonon scattering, for the density parameter 
ao =0.1 and three temperatures t p = 0.1, 1, 10 in bcc, 
fee and hep lattices. The curves for fee and hep crystals 
almost coincide. Dots show the effective structure factor 
which corresponds to the one-phonon approximation 



Figure g shows the dependence of S c s (q) on y at ao = 
0.1 and three temperatures t p = 0.1, 1, and 10 for three 
types of Coulomb crystals, bcc, fee, and hep. The results 
are seen to be almost independent of lattice type, and 
the curves for fee and hep lattices are indistinguishable. 
If tp 1, the one-phonon approximation appears to be 
highly inaccurate for y ^ 0.4. 

Note that the present treatment of phonon scattering 
is valid as long as T ^ T u , where T u ~ T p Z 1 / 3 e 2 / (Hv F ) 
is the temperature below which umklapp processes are 

c 



frozen out (see, e.g., Raikh & Yakovlev 1982), and vf 



Ses(q) = (• 



1 e 



-2W 



(28) 



is the electron Fermi velocity. If T ^ T u the electron- 
phonon scattering can be described in the free electron ap- 
proximation, and the main contribution to the scattering 
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comes from umklapp processes (Yakovlev & Urpin 1980 



Raikh & Yakovlev 1982). The density dependence of T u is 
shown in Fig. 1 of the article by Baiko & Yakovlev ( 1996[ ); 
T u is so low that a study of neutrino emission at lower 
temperatures is of no practical importance. 

3.3. Static-lattice contribution 

This contribution corresponds to neutrino emission due to 
Bragg diffraction of electrons in a crystal. It was widely 
assumed for a long time that the process could be stud- 
ied in the free-electron approximation (e.g., Flowers 1973; 
Itoh ct al. 1984a). However, recently Pethick & Thors- 



son (1994) have pointed out the importance of electron 
band-structure effects. The gaps in the electron disper- 
sion relation at the boundaries of the Brillouin zones can 
reduce noticeably the static-lattice contribution as com- 
pared to that obtained in the free-electron approximation. 
The general formalism for evaluating the static-lattice 
contribution for strongly degenerate relativistic electrons 
with proper treatment of band-structure effects has been 
developed by Pethick & Thorsson (|1997[). Their Eq. (28) 



is valid for spherical and nonspherical nuclei and can be 
written as 



27rG|fc F C; 
567n 9 c 8 



■{ksTfj 



1.254 x 10 9 (pi 2 Y e ) 1/3 7f J erg s" 1 cnT 



where 



J = w I (y> t v) 



(30) 



(31) 



is a dimensionless function given by a sum over all re- 
ciprocal lattice vectors K ^ for which K/2 lies within 
the electron Fermi surface; y = K/(2kp) (with y < 1). 
The function I(y,ty) is given by Eq. (29) in Pethick & 
Thorsson (1997), which can be rewritten as 



i{y,tv) = — 

8tt 7 y\ 



dzi dz2 



x[(w 2 - w)(w - wi)] 3/2 , 



dw- 



(32) 



where y± = yl 



E,.> = \/l + (*iTy) • 



E 1 +E 2 Ty±V(E 1 +E 2 y-z^ 
t v y 2 



(33) 



ty is defined in terms of the Fourier transform of the lat- 
tice potential (0) with q=K by 

J_ = \Vk\v± = ^e Pz \F(K)\ W[K) 

t v k B T K 2 k B T \e(K)\ y±G ' 10 1 

and other notations were introduced in Eq. (Q). 



3r 2 
1 

12Z 



(35) 



For a lattice of spherical nuclei, we can use Eq. (13) 
with 

i.L = ^(9»Z/4)" 3 J 

E 4 !W ..^w. 

y H K )\ 

According to Eqs. (|), @ and @, the Debye- Waller 
factor suppresses the electron-lattice interaction at large 
reciprocal lattice vectors K and weakens the neutrino 
emission. Computing Q s i directly as a sum of 3D integrals 
( |32"| ) is time consuming since the number of reciprocal lat- 
tice vectors K involved is generally large (~ AZ terms). 
We simplify computation by producing an analytic fit to 

i{y,t v ). 

Analytical asymptotes of I(y, ty) can be derived in the 
limiting cases of high and low temperatures. In the high- 
temperature limit, ty 3> 1, Pethick & Thorsson ( |1994[ , 
1997D obtained 

^-(l + ^f-lnyV (36) 



I 



ylv 



Inserting this asymptote into Eqs. (Ba) and (O) one immc 



diately reproduces the well-known result of Flowers (1973) 
and Itoh et al. (1984a) for the static-lattice contribution 
when band-structure effects are neglected. Replacing the 
sum over K by an integral over q, we arrive at the expres- 
sion 

nl , \F(q)\ 2 e- 2W ( 2y 2 
dy 

'yo 



-(o) 



1 



lny 



(37) 



which is similar to Eqs. ([l6]) and ( pl| ) in the liquid and 
for the phonon contribution in the solid, respectively. The 
Debye-Waller exponent e~ 2W is seen to play the role of 
the diffraction part of an "effective static structure factor" 
that defines the static-lattice contribution (smoothed over 
diffraction peaks due to replacing summation by integra- 
tion). Thus the sum Lp^ + L^ in a crystal can be written 
in the same form ( |l6| ) as Lii q in a liquid, with an effective 
structure factor S so i(q) — e~ 2W + S c s(q). We have veri- 
fied that S so i(q) resembles the structure factor S(q) in a 
strongly coupled liquid (Young et al. 1991) for ion coupling- 
parameters 100 r ^ 225 if we smear out the familiar 
diffraction peaks in S(q); the integral contributions of the 
two factors are nearly the same. This elucidates the sim- 
ilarity of neutrino-pair bremsstrahlung in a liquid and a 
crystal (Sect. 4). 

The asymptote ( |3^ ) is temperature-independent. The 
lowest-order thermal correction can be taken into account 
by introducing the factor [1 + (63/40)(ntyy)~ 2 ]. 

The low-temperature asymptote for ty -C y±/y 2 is 



(Pethick & Thorsson |1997| ) 



i{y,t v ) = 



189 



27T n /2(l 



x exp 



y±) 1/2 (i- 

2 

'tv0- + y±) 



yx) 2 (yxtyf/ 2 



(38) 
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while for yi/|/ 2 < C 1 the latter asymptote has to be 
multiplied by y/n[y x / \t V y 2 )f /2 ■ 

We calculated I(y,ty) from Eq. (|3^) numerically for 
wide ranges of y and ty and derived an analytic fit which 
reproduces numerical results and the asymptotes: 



i(y,tv) 



0.3088 (1 + 8A16y±t v ) 
y(l + yx) 3 ' 2 (y±t v )*/*D 
1 ( ^ 0.4031 



exp 



VxD 



1 



t v y 2 + OMv + 0.2678 



1 



2y\ 
— In y 

V± 



exp 



t v H 



(39) 



where D = u 
1.772, H = 1 



y 



0.7124m 4 - 1.689m 3 + 5.237 m 2 -0. 2u 

212 tyy 2 , and u = yyjty /y±. The 



mean error of the fit is 2%, and the maximum error of 
7.4% occurs at t v = 0.045 and y = 0.987. The fit permits 
a rapid and accurate evaluation of the static-lattice con- 
tribution from Eqs. ([50]) and (|3l]) or ( |l3"| ) and (|3^) for any 
plasma parameters of practical interest. 

For a lattice of spherical nuclei, we can use Eq. ( |35| ) 
with the Debye-Waller factor and the nuclear form fac- 
tor, and calculate the sum over reciprocal lattice vectors 
for bcc, fee and hep crystals. In the case of non-spherical 
nuclei, we use the more general Eqs. ( pO] ) and (|l|) includ- 
ing the form factor but setting W — 0, since the Debye- 
Waller factor is unknown. The sum over reciprocal lattice 
vectors in Eq. ( |3l| ) for a lattice of non-spherical nuclei 
is different from that for ordinary crystals. In the case of 
rod-like phases, the lattice is the simplest two-dimensional 
triangular one, for slab- like nuclei the lattice is purely one- 
dimensional, and in the case of the Swiss-cheese phase with 
neutron drops immersed in nuclear matter we shall assume 
that the lattice is bcc. 

4. Results and discussion 

Let us outline the main properties of neutrino-pair 
bremsstrahlung by relativistic degenerate electrons in 
neutron-star crusts. The results of Sect. 3 allow us to cal- 
culate the neutrino emission for bcc, fee and hep Coulomb 
crystals of atomic nuclei. However, in all three cases the 
emission appears to be practically the same because of the 
close similarity of the crystals. For simplicity we therefore 
consider bcc crystals throughout this section. 

Figures |^ and ^ show the temperature dependence 
of the normalized neutrino emissivity L for two densi- 
ties, p = 10 10 and 10 12 g cm~ 3 , in the outer and in- 
ner crusts, respectively. Figure [3] is plotted for matter 
composed of iron (Z = 26, A = 56), while Fig. H cor- 
responds to the ground-state matter composed of 49° Zr 
nuclei (Z = 40, A = 250, where A is the number of nucle- 
ons per Wigner-Seitz cell, see Negele & Vautherin 1973| ). 
Vertical dotted lines separate liquid and solid phases. The 
upper (dashed) line in the liquid phase is obtained from 
Eq. (gj) with the non-Born corrections included in the 




lg T [K] 



Fig. 3. Temperature dependence of the normalized neu- 
trino emissivity L for iron matter at p = 10 10 g cm~ 3 . 
Solid lines: Born results for the liquid phase; phonon 
and static-lattice contributions, as well as the total func- 
tion ( fl4| ) for the crystalline phase. Dotted line: the to- 
tal function for crystalline phase but without band- 
structure effects. Dashed line: non-Born result in the liq- 
uid phase. Dot-dashed line: one-phonon approxi matio n for 
the phonon contribution (Yakovlev & Kaminker 1996 ). All 
curves but one in the liquid phase are obtained in the Born 
approximation 



function i?NB (Haensel et al., 1996j ). The lower (solid) line 
is also obtained from Eq. ( [l5| ) but in the Born approxima- 
tion (i?NB = !)■ Solid lines in the crystalline phase show 
the phonon contribution (Sect. 3.2), the static-lattice con- 
tribution (Sect. 3.3), and the total function L so i given by 
Eq. ([lj]). The dotted line also gives the total function but 
neglecting the band structure effects in the static-lattice 
contribution (by using the asymptote ( |36| ) in Eq. (|35|)). 
Finally, the dot-dashed line displays the phonon contribu- 
tion calculated in the one-phonon approximation adopted 
in previous articles (Sect. 3.2). 

The temperature profiles of L in Figs. || and ^ are sim- 
ilar. The phonon contribution is generally several times 
smaller than the static-lattice one. Each term in the static- 
lattice sum ([}5]) is suppressed exponentially with decreas- 
ing temperature but the sum itself decreases more like 
a power law because, for the smallest reciprocal lattice 
vectors \K\, the exponential decrease of the contribution 
starts to operate at much higher temperature than that 
for la rger reciprocal lattice vectors (Pethick & Thorsson 
1997 ). At very low temperatures, the reduction of con- 
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Fig. 4. Same as in Fig. |^ but for ground-state matter 
composed of 2 4oZr nuclei at p = 10 12 g cm -3 



tributions from all reciprocal lattice vectors becomes ex- 
ponential, and the total static-lattice contribution is sup- 
pressed exponentially. However, such low temperatures are 
of no practical importance. Generally, the static-lattice 
neutrino emission is usually partially suppressed by band- 
structure effects, and these effects are substantial. 

All calculations of neutrino-pair emission in the crys- 
talline phase have been made in the Born approximation. 
On the other hand, Haensel et al. ( 1996 ) calculated L\i q 
beyond the Born approximation. For comparison of the re- 
sults in the crystalline and liquid phases, we present also 
the Coulomb logarithm Li; q determined in the Born ap- 
proximation. This curve matches the function L so i in the 
solid phase much better and makes L an almost continuous 
function of temperature at the melting point. We believe 
that the non-Born curve in the solid phase (which is dif- 
ficult to calculate exactly) would match equally well the 
non-Born curve in the liquid phase, because of the similar- 
ity of the bremsstrahlung in a liquid and in a crystal, men- 
tioned in Sect. 3.3. The state of a Coulomb system (liq- 
uid or solid) is expected to have little effect on neutrino- 
brcmsstrahlung, since in both cases neutrino emission oc- 
curs due to scattering of electrons by fluctuations of elec- 
tric charge produced by ions (nuclei). In the solid phase 
near the melting point both the phonon and the static- 
lattice contributions must be included for the total emis- 
sion rate to match the results in the liquid phase. The one- 
phonon approximation is seen to be generally quite accu- 
rate at low temperatures, T <C T m (actually at T <§; T p , as 
discussed in Sect. 3.2) but underestimates the phonon con- 



tribution near the melting point. It is the proper inclusion 
of multi-phonon processes that makes the phonon contri- 
bution larger and almost removes the jump of the total 
neutrino emissivity at the melting point. Notice that the 
phonon and static-lattice contributions at the melting be- 
come nearly equal. The physical properties of a Coulomb 
liquid and a Coulomb solid near the melting point are 
nearly the same. To verify this statement we have shifted 
artificially the melting temperature (which actually cor- 
responds to r = r m = 172, Sect. 2). We have considered 
the cases of a supercooled liquid by taking r m = 225, 
f or wh ich the liquid-state structure factors of Young et al. 
( 1991 ) are available, and the case of a superheated crystal 
by taking r m = 100. In all the cases the total normalized 
neutrino emissivities L do not differ noticeably from those 
obtained at r m = 172, and the discontinuity of L at the 
melting point is minor. 

Figure || compares the present results for 56 Fe matter 
at p = 10 9 g cm -3 with the familiar results by N. Itoh and 



his group reviewed recently by Itoh et al. ( 1996Q and with 
the results by Yakovlev & Kaminker (1996) for the phonon 
contribution obtained in the one-phonon approximation^] 
The temperatures displayed are rather low, so the one- 
phonon approximation almost coincides with the multi- 
phonon calculation. One can see a transition from power- 
law to exponential decrease in our static-lattice curve with 
decreasing temperature at T ~ 10 7 K. For lower T, the 
phonon contribution dominates over the static-lattice one. 



The static-lattice contribution given by Itoh et al. (1996) 



is underestimated by several orders of magnitude. The au- 
thors calculated this contribution neglecting the electron 
band structure and multiplied this result by a factor which 
should approximately take into account the suppression of 
the neutrino emission due to band structure. The latter 
factor was chosen on the assumption that one particular 
reciprocal lattice vector gave the dominant contribution 
at all temperatures. If this were the case, the suppres- 
sion would be exponential. However, as shown above, the 
suppression of contributions from the various reciprocal 
lattice vectors sets in at different temperatures, and the 
resulting suppression of the total rate is much weaker than 
given by the approximation of Itoh et al. (1996). In fact, 
the approximate suppression factor introduced by Itoh 



et al. (1996) makes the static-lattice contribution negli- 
gible at all densities and temperatures of practical inter- 
est. Since the actual static-lattice emission is commonly 
several times larger than the phonon one, the approach 



1 Note that Yakovlev & Kaminker (1996), while comparing 
their results for the one-phonon contribution in the iron and 
carbon plasmas at p = 1 9 g c m" 3 (their Figs. 3 and 4) with 



the results of Itoh et al. (1996), inaccurat ely plo tted (by long 
dashes) the fit expressions of Itoh et al. (1996). The c orrect 
curves are closer to the results by Yakovlev & Kaminker ( 1996 ) 
and are plotted in the present Figs. ^ and []. Two of the authors 
(DY and AK) are grateful to Prof. N. Itoh for pointing out this 
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Fig. 5. Phonon (ph), static-lattice (si) and total normal- 
ized neutrino emissivities £ p h, L s \, L so \ vs T (lower hori- 
zontal scale) or t p (upper horizontal scale) for a crystal of 
56 Fe nuclei at p = 10 9 g cm -3 calculated in the Born ap- 
proximation. Solid curves show the present results, dashes 
show the phonon contribution i n the one-phonon approx- 
imation (Yakovlev & Kaminker 
results of Itoh et al. 



0L996D 



1996) while dots are the 



of Itoh et al. ( 1996| ) systematically underestimates the to- 
tal neutrino emissivity in the solid phase of dense matter. 
In addition, the phonon contribution given by Itoh et al. 
( |1996| ) appears to be overestimated for T ^ 10 7 K ( Fig, ph 
as a result of the analytic fits proposed by Itoh et al. (1996) 
being insufficiently accurate (Yakovlev & Kaminker 199(f ). 

Figure ^ presents a similar comparison of the results 
but for carbon. The case of carbon at high density is 
extreme since zero-point vibrations of the light carbon 
ions become very strong. Nuclear reactions and beta cap- 
tures tend to transform carbon into heavier elements. The 
Debye- Waller factor is very large owing to zero-point vi- 
brations (a oc u-i/y/A{Z at t p <C 1, see Eq. (||)) 



suppresses drastically the static-lattice contribution and 
makes it generally smaller than the phonon contribution. 
The simplified treatment of the band-structure effects by 



Itoh et al. (1996) damps the static-lattice contribution es- 
pecially strongly, by several orders of magnitude, making 

Figure shows the density dependence (10 9 g cm~ 3 < 
p < 10 13 g cm~ 3 ) of the neutrino emissivity at three val- 
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Fig. 6. Same as in Fig. | but for carbon crystal 



ues of T for the ground-state and accreted matter. Here 
and in all subsequent figures the emissivities are calcu- 
lated in the Born approximation and multiplied by the 
non-Born correction factor ( p"8| ) as discussed in Sect. 3.1. 
Note that self-consist ent m odels of accreted matter (e.g., 
Miralda-Escude et al. |l990|) correspond to T ~ 10 8 K. We 
use the accreted model for higher T to illustrate how vari- 
ations of nuclear composition affect the neutrino emission. 
For a particular temperature and density, matter will be 
either liquid or solid. In the liquid state, we calculate Q 
using Eq. (|l5|). In the solid state, we present the total 
and phonon neutrino emissivities for ground-state mat- 
ter and the total emissivity for the accreted matter. In 
accreted matter the ratio of the total emissivity to the 
phonon one is qualitatively the same as in ground-state 
matter. We show also the total neutrino emissivity for 
ground-state matter neglecting band-structure effects. In 
the displayed density range, matter is entirely solid for 
T = 10 8 K; there is one melting point for T = 10 8 - 5 K 
It (lg/5 m [gem -3 ] = 9.17, for ground-state matter) which 



separates liquid (at p < p m ) and solid (at p > /? m ); 
and there is a series of melting points at T = 10 9 K 
due to the non-monotonic behaviour of the melting curves 
T m = T m (p) associated with strong variations of the nu- 
clear composition. The positions of the melting points can 
be deduced from Fig. [I]. For ground-state matter, these 
positions can also be traced in Fig. Q from the appear- 
ance of the phonon contribution (dashed lines). With de- 
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Fig. 7. Density dependence of the neutrino 
bremsstrahlung emissivity at T = 10 8 , 10 8 5 and 
10 9 K for ground-state and accreted matter calculated 
with the form factor (|ll|) . Solid and dashed lines show the 
total and phonon emissivities, respectively, for ground- 
state matter; dots present the same total emissivity but 
obtained neglecting the electron band-structure effects. 
Dot-dashed lines display the total emissivity for accreted 
matter 



Fig. 8. Density dependence of the neutrino 
bremsstrahlung emissivity at T = 3 x 10 8 K for 
ground-state matter calculated with the form factor (pi]). 
Solid line: total emissivity. Also shown are: static-lattice 
contribution (si), phonon contribution (ph) and the total 
contributions obtained neglecting either band-structure 
effects (no band) or Debye-Waller factor (no DW) or 
nuclear form factor (no ff) 



creasing temperature the solidification front shifts toward 
lower densities. The accreted matter solidifies at higher 
densities than does ground-state matter due to the lower 
atomic number. At the melting points the neutrino emis- 
sivities exhibit some jumps, but these are small because, 
as we remarked earlier, neutrino bremsstrahlung does not 
change qualitatively while passing from liquid to solid 
matter. Other, stronger jumps are associated with varia- 
tions of the nuclear composition (Sect. 2). Notice that the 
jumps of both types may be ignored in practical applica- 
tions. The reduction of the neutrino emission by the band- 
structure effects becomes stronger with decreasing tem- 
perature and reaches one order of magnitude for T ~ 10 8 
K and p <; 10 11 g cm~ 3 . The band-structure reduction is 
power-law (non-exponential) for the parameters displayed 
in Fig. 0. The ratio of the phonon contribution to the 
static-lattice one remains nearly constant for a wide range 
of temperatures much below the melting temperature, and 
the static-lattice contribution is several times larger than 
the phonon one. The neutrino bremsstrahlung in the ac- 
creted matter is lower than in the ground-state matter 
due to the lower atomic number, but the difference is not 
large. 



The effects of various physical factors on the density 
dependence of the bremsstrahlung emissivity in the lattice 
of ground-state matter at T = 3 x 10 8 K is shown in 
Fig. H We present the total emissivity, and also the static- 
lattice and phonon emissivities. In addition, we show the 
total emissivity calculated neglecting either the effects of 
band-structure, the Debye-Waller factor, or the nuclear 
form factor. Here by the phrase "neglecting the Debye- 
Waller factor" we mean the one-phonon approximation in 
which the Debye-Waller factor is set to zero. One can see 
that the phonon contribution is noticeably smaller than 
the static-lattice one, and the ratio of the static-lattice 
and phonon contributions increases slowly with density, 
reaching a factor of about 5 at p — 10 13 g cm~ 3 . The 
effect of the form factor also increases with density while 
the effect of the Debye-Waller factor becomes lower. 

The neutrino emissivities presented in Figs. |3|-|q are 
calculated with the simplified atomic form factor ( |11[ ) ap- 
propriate for a step-like proton distribution within the nu- 
clei. We have verified that the simplified form factor gives 
quite accurate results for p <, 10 13 g cm -3 . However, at 
higher p the neutrino emissivity becomes sensitive to the 
shape of the proton distribution. Then the form factor 
based on the proton-density distribution (|l^) seems to be 
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Fig. 9. Density dependence of the bremsstrahlung emis- 
sivity at T = 3 x 10 8 K for the ground-state matter com- 
posed of spherical nuclei. Long-dashed lines are calculated 
with the step-like form factor (ff ) given by Eq. (|lTJ) for the 
nuclear composition by Negele & Vautherin (NV,1973). 
Short-dashed lines are calculated with the same form fac- 
tor but for the smooth-composition (SC) model of matter 
(Sect. 2). Solid lines correspond to the SC model and more 
realistic form factor for the proton distribution (|l2]). Up- 
per lines show the total emissivities, while the lower ones 
give the static lattice contribution alone 



more reliable, as discussed in Sect. 2. We have made a 
series of calculations with this more realistic form factor 
making use of the smooth-composition (SC) model of the 
ground-state matter (Sect. 2). Figure |9| is an extension of 
Fig. H to higher densities (10 12 g cm~ 3 < p < 8 x 10 13 g 
3 ) for matter containing spherical nuclei. The highest 



cm 

density, p = 8 x 10 i3 g cm -3 , is close to the transition 
to the phases with nonspherical nuclei. We now compare 
the neutrino emissivities calculated with the realistic and 
step-like form factors. The realistic form factor is included 
in both the phonon contribution (through Eq. d2l|)) and 
the static-lattice contribution (Sect. 3.3). The emissivities 
obtained with the realistic form factor for the SC model 
are compared with those obtained with the step-like form 
factor for the SC and Negele- Vautherin models of mat- 
ter. The total emissivities (upper curves) are seen to be 
close to the static-lattice emissivities (lower curves) , which 
indicates that the static-lattice contribution is dominant. 
The emissivities in the Negele- Vautherin model show the 
familiar jumps (e.g., Figs. ^ and |^) associated with vari- 
ations of the nuclear composition. After averaging over 



these jumps, they reproduce the emissivities derived in 
the SC model (with the same form factor). The emissivi- 
ties obtained with the realistic and step-like form factors 
are seen to be close for p ^ 10 13 g cm~ 3 . However, when 
the density becomes higher than 10 13 g cm~ 3 the realistic 
form factor decreases the electron-nucleus interaction and 
reduces noticeably the neutrino emission compared with 
the case of the step-like form factor. 
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Fig. 10. Density dependence of the bremsstrahlung emis- 
sivity at T = 3 x 10 8 K for five phases of spherical and 
nonspherical atomic nuclei at the bottom of a neutron star 
crust for the SC model of matter calculated with the real- 
istic form factor. Long-dash line is obtained by assuming 
the nuclei to be spherical to the crust bottom. Short-dash 
line is calculated for all phases neglecting the band struc- 
ture, while dotted line is obtained neglecting the band 
structure and the form factor 

Figure |l^ is an extension of Figs. - 1| towards higher 
densities. It shows the density dependence of the neu- 
trino bremsstrahlung emissivities through all five phases 
of spherical and nonspherical nuclei for T = 3 x 10 s K. 
In order to display all the phases we use a linear density 
scale rather than a logarithmic one; pu is the density in 
units of 10 14 g cm" 3 . Densities in excess of 1.43 x 10 14 g 
cm~ 3 correspond to uniform matter in the neutron star 
core. The emissivities are calculated using the SC model. 
In the spherical phase, the phonon contribution is included 
and the Debye- Waller factor is taken into account. In the 
nonspherical phases, we neglect the Debye-Waller factor 
and the phonon contribution (see above). This circum- 
stance is partly responsible for the jumps in the neutrino 
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emissivities at p — 9.73 x 10 13 g cm -3 , the interface be- 
tween the phases with spherical and cylindrical nuclei. All 
the curves are calculated with the realistic form factor. 
The long-dashed curve is obtained assuming the nuclei 
remain spherical up to the highest densities at the very 
bottom of the neutron star crust and by including the 
phonon contribution and the Debye-Waller factor. The 
parameters of such nuclei are appropriate to the Negele- 
Vautherin model of dense matter, smoothed over jumps. 
We also present the emissivities calculated for all phases 
of nonspherical nuclei neglecting either the band-structure 
effects, or both the band structure and the form factor. 
The neutrino emission at p ~ 10 g cm -3 is very sen- 
sitive to the proton charge distribution. A neglect of the 
form factor leads to overestimation of the neutrino emis- 
sivity by 1 - 1.5 orders of magnitude. The effects of non- 
spherical phases are also rather important. Non-sphericity 
of the nuclei mainly lowers the neutrino emission by re- 
ducing the dimension of the sums over reciprocal lattice 
vectors in Eq. ([3l]). The reduction can exceed one order 
of magnitude (cf the solid and dashed curves near the 
cylinder-slab interface). More work is required to calcu- 
late the Debye-Waller factor and the phonon contribution 
and determine accurately the bremsstrahlung emission for 
nonspherical nuclei. However the spherical-nucleus model 
probably represents a reasonable upper limit to this emis- 
sivity for densities where the nuclei can be nonspherical. 

Figure [ll] shows the density dependence of neu- 
trino bremsstrahlung for the ground-state matter of the 
neutron-star envelopes at six temperatures, from 10 8 K 
to 1.8 x 10 9 K, in the density range 10 10 g cm " 3 < p < 
1.4 x 10 14 g cm -3 . The p — T domain displayed is the most 
important one for application to neutron-star cooling. The 
dotted curves are calculated using the Negele-Vautherin 
model of matter and the step- like form factors (fTl|). The 
dashed lines are obtained for the SC model with the re- 
alistic form factor and on the assumption that nuclei are 
spherical to the bottom of the crust. The solid lines are 
also derived for the SC model with the realistic form factor 
but with allowance for the phases of nonspherical nuclei. 
Sharp decreases of the curves at p 10 14 g cm" 3 occur 
because the nuclei become nonspherical (see Fig. [jl]) . 

Our calculations for spherical nuclei with the realistic 
form factor in wide density and temperature ranges, 10 9 g 
cm 

K can be fitted by the expression 



3 < P < 1.5xl0 14 gcm" 3 and 5 x 10 7 K < T < 2xl0 9 



lgQ [erg cm" 3 s" 1 ] = 11.204 + 7.304 r + 0.2976 r 

-0.370 t 2 + 0.188 tt - 0.103 r 2 + 0.0547 r 2 r 
-6.771g(l + 0.228p/p ), (40) 

where r = lgT 8 , r = \gpii, and p = 2.8 x 10 14 g cm" 3 
is the standard nuclear-matter density. The relative error 
of this fit formula generally does not exceed 1% (in lg Q) 
over the indicated p — T domain. The accuracy of the fit is 
seen from Fig. [ll]. The fitting formula reproduces the main 
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Fig. 11. Density dependence of the bremsstrahlung emis- 
sivity for the ground-state matter of the neutron-star crust 
at six temperatures T in the model by Negele & Vautherin 
(NV, dots) Q1973 ) with the step-like form factor (ff), and 
in the SC model with the realistic form factor assuming 
either the nuclei are spherical to the crust bottom (dashes) 
or including non-spherical phases (solid lines). Filled cir- 
cles show our fits (^0|) to dashed lines 

features of the neutrino bremsstrahlung radiation in a very 
wide density domain 10 9 g cm" 3 < p < 10 14 g cm" 3 , 
where the atomic nuclei are expected to be spherical, and 
it probably gives a realistic estimate of the emissivity for 
higher densities, 10 14 g cm" 3 < p < 1.4 X 10 14 g cm" 3 , 
where the nuclei become nonspherical. This fit formula 
should be quite sufficient for many applications. 

Finally, in Figs, [l^ and [ll^ we compare neutrino 
bremsstrahlung by electrons in the neutron-star crust with 
other neutrino emission mechanisms for T = 10 9 K and 
10 8 K, respectively. The other mechanisms considered are 
neutrino emission due to plasmon decay (e.g., Itoh et al. 



1996) and due to the nucleon Cooper-pair formation un- 
der the action of nucleon superfluidity (e.g., Yakovlev et 
al. 1998). The SC model of ground-state matter is used. 



In principle, some contribution at low T may come from 
neutrino emission due to scattering of electrons by charged 
impurities (e.g., Haensel et al. 1996] ). We ignore this mech- 
anism here since it is determined by impurity parameters 
which are largely unknown. We show also the neutrino 
emission produced by some mechanisms in uniform mat- 
ter of the neutron-star core: the standard neutrino en- 
ergy losses and the nucleon Cooper-pair neutrinos in su- 
perfluid uniform matter. It is assumed that uniform mat- 
ter has a moderately stiff equation of state proposed by 
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Fig. 12. Density dependence of the neutrino emissivi- 
ties produced in a neutron star crust and core at T = 
10 9 K by various neutrino generation mechanisms: elec- 
tron bremsstrahlung (e-brems), plasmon decay (plasma), 
Cooper pairing of nucleons (Cooper pairs) in the mod- 
els of strong (s) and weak (w) neutron superfluidity (SF), 
and standard neutrino generation mechanisms (N-stand.) 
in strongly (s) and weakly (w) superfluid uniform matter 
(see text for details) 



Prakash et al. ( 1988) ( the same version as was used by 
Page k Applegate F992Q . The standard neutrino processes 
include the modified Urea reactions and nucleon-nucleon 
bremsstrahlung. All the standard reactions are partially 
suppressed by the combined action of the neutron and 
proton supcrfluids in uniform matter. The reaction rates 
and suppression factors are taken in the form described 
by Levenfish & Yakovlev (1996). The proton and neutron 



superfluid critical temperatures depend on density. We as- 
sume singlet-state pairing of protons, and either singlet- 
state or triplet-state pairing of neutrons. Neutron pairing 
is expected to occur in a singlet state in matter at densities 
less than roughly that of the core-crust interface, and in a 
triplet state at higher densities. For any density we adopt 
the type of neutron pairing that corresponds to the higher 
critical temperature. The strong density dependence of the 
neutrino emissivities in uniform matter (Figs. |l^ and |l3| ) 
is due to the pronounced density dependence of the neu- 
tron and proton critical temperatures at p ~ po- Since 
critical temperatures are very sensitive to the microscopic 
model of the nucleon interaction, we have considered two 
cases, corresponding to strong (s) and weak (w) nucleon 
superfluids (SFs). The strong superfluid model is based on 
the rather large superfluid gaps calculated by Elgar0y et 



Fig. 13. Same as in Fig. |T| but for T = 10 8 K. Plasmon 
decay, standard neutrino emission from uniform mater and 
Cooper-pairing emission in the strongly superfluid matter 
become negligible 



al. (1996) for singlet-state pairing (with maximum gap of 



about 2.5 MeV as a function of nucleon number density) 



and by Hoffberg et al. (1970) for triplet-state pairing. The 



weak superfluid model makes us e of t he small superfluid 
gaps derived by Wambach et al. ( 1993 ) (with a maximum 
gap of about 1 MeV) for th e singlet superfluid and by 
Amundsen & 0stgaard ( |1985| ) for the triplet neutron su- 
perfluid. For singlet pairing we regard the weak pairing 
case to be the more realistic because the calculations of 



Wambach et al. (1993) included the effects of induced in- 
teractions. 

The mechanism of neutrino production due to Cooper 



pairing of nucleons was proposed by Flowers et al. ( 1976|) 
and independently by Voskresensky & Senatorov (19861 



1987) for uniform matter. A critical comparison of these 



works has been done by Yakovlev et al. (1998) who con- 
sidered also the case of triplet neutron pairing. The the- 
ory predicts a powerful maximum of the Cooper-pair neu- 
trino emission when the temperature falls below the crit- 
ical temperature T c for neutrons or protons. However, at 
still lower temperatures, T -C T c , the emission falls ex- 
ponentially The emission is much stronger for neutrons 
than for protons due to the smallness of axial-vector elec- 



troweak currents for protons (Yakovlev et al. 1998). In 
spite of its long history, this process was forgotten for a 
long time, and it has been included in cooling simulations 



only rece ntly (P age |1997| , |1998[ Scha ab et al. |1997| ; Leven- 
fish et al. |1998| ; Yakovlev et al. |l998| ). We make use of the 
results by Yakovlev et al. ( 1998 ) for uniform matter and 
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for the crust. In the uniform matter, we take into account 
neutrino emission due to pairing of neutrons and protons. 
In the crust, the Cooper neutrino emission is evaluated 
including the contribution from free nucleons alone (from 
free neutrons in all phases of matter in the inner crust, and 
from free protons in the phases with cylinders and spheres 
of neutrons surrounded by nuclear matter). In principle, 
there can be a substantial contribution (Yakovlev et al. 
199S) from the non-uniform distribution of the nucleons 
within atomic nuclei (Sect. 2). We neglect this effect here 
but we intend to discuss it in a separate article. 

In a high temperature plasma, T = 10 9 K, at p 10 11 
g cm -3 (Fig. [l2|) , the process most competitive with elec- 
tron bremsstrahlung is the well-known plasmon decay 
(Itoh et al. |l996 ). However its rate falls exponentially 
with decreasing T, and the process almost dies out at 
T = 10 8 K (Fig. |i~3| ). The standard neutrino emission 
from uniform matter is greatly reduced by the nuclcon 
superfluidity. It decreases exponentially when the tem- 
perature becomes much smaller than the superfluid crit- 
ical temperatures of nucleons. For instance, the standard 
emission is quite substantial at T = 10 9 K but becomes 
much less significant at T = 10 s K (practically negligi- 
ble for the case of strong superfluid). Neutrino emission 
by Cooper pairing in the uniform matter and in the in- 
ner crust is also exponentially suppressed when the tem- 
perature is much lower than the critical temperatures of 
neutrons and protons. Accordingly the temperature and 
density dependence of the Cooper-pair neutrino emissivity 
is very strong. If T = 10 9 K and the superfluidity is strong, 
we have two peaks of Cooper-pair neutrinos: one near the 
neutron drip point (at p ~ 10 12 g cm -3 ), and a very nar- 
row peak near the core-crust interface (p « 1.4 x 10 14 g 
cm -3 ). Both peaks are produced by neutron pairing, and 
are pronounced since the neutron critical temperature is 
sufficiently small (only slightly exceeds T) in the indicated 
density ranges even for strong superfluidity. The small- 
ness of T c in the first density range is associated with a 
low number density of free neutrons near the neutron drip 
point, and the smallness in the second density range cor- 
responds to the transition from singlet to triplet neutron 
pairing. Outside these density ranges the neutron critical 
temperature is too high and the emissivity of the Cooper 
neutrinos is exponentially suppressed. When the temper- 
ature decreases the emissivity becomes smaller, and the 
process dies out at T = 10 8 K (cf. Figs. || and |l3|). 

If T = 10 9 K and the superfluidity is weak, the Cooper 
pairing appears to be the dominant neutrino emission pro- 
cess in a large fraction of the neutron star crust since the 
neutron critical temperature is not much higher than T. 
However, the process is practically switched off at low den- 
sities p < 2 x 10 12 K, because a weak neutron superfluid 
has not yet occurred at these densities (T c < T = 10 9 
K). When the temperature drops to 10 8 K, the neutrino 
emission due to Cooper pairing is suppressed. Neverthe- 
less, two high peaks of the emissivity survive (similar to 



those for a strong superfluid at T = 10 9 K). The first one 
corresponds to low p ~ 10 12 g cm" 3 (T c is not much higher 
than 10 8 K), and the second corresponds to p ~ 1.4 x 10 14 
g cm -3 , where T c is low, which corresponds to the transi- 
tion from a singlet to a triplet neutron superfluid. 

We conclude that the main contribution to neutrino 
emission from a neutron-star crust comes from two pro- 
cesses, the neutrino-pair bremsstrahlung and Cooper pair- 
ing of neutrons. The bremsstrahlung neutrino emission has 
been calculated rather reliably, excluding possibly in the 
phases of nonsphcrical nuclei near the core-crust interface. 
The mechanism operates in a wide ranges of densities and 
temperatures, and the density dependence of the emis- 
sivity is generally smooth. The neutrino emission due to 
the Cooper pairing of nucleons is extremely sensitive to 
the model adopted for the superfluid gaps in the nucleon 
spectra. This mechanism is more important for lower gaps 
(weaker superfluid). The emissivity can be a sharp func- 
tion of density and temperature. We remark that the mi- 
croscopic models that correspond to weak superfluidity are 
likely the more reliable since they incorporate the effects 
of induced processes (screening) in the effective neutron- 
neutron interaction. In addition, in the phases with non- 
spherical nuclei, the neutron superfluid gap in the matter 
is expected to be reduced by the presence of nuclei, in 
which the matter has a higher density, and a smaller pair- 
ing interaction, than in the neutrons outside nuclei. We 
expect to consider the Cooper-pair process in more detail 
in a future publication. 

5. Conclusions 

We have analysed the neutrino pair emission (|]J) due to 
the bremsstrahlung of degenerate relativistic electrons at 
densities from 10 9 g cm -3 to 1.5 x 10 14 g cm -3 and tem- 
peratures from 5 x 10 7 to 5 x 10 9 K in the neutron-star 
crusts. We have presented the expressions for the neu- 
trino emissivity from a plasma of liquid and solid atomic 
nuclei taking into account the effects of finite sizes (the nu- 
clear form factor) of the nuclei. In solid matter, we have 
studied the static-lattice and phonon contributions to the 
neutrino bremsstrahlung with allowance for the electron 
band-structure and multi-phonon processes. We have con- 
sidered bcc, fee and hep Coulomb crystals, and showed 
that the neutrino emission is insensitive to the lattice type. 
We have made use of two models of matter in the neu- 
tron star crusts: ground-state matter and accreted matter. 
We have proposed a smooth-composition model of ground- 
state matter to analyse the neutrino bremsstrahlung near 
the bottom of the crust, where the shapes of the local 
nuclcon density distributions over Wigner-Seitz cells be- 
come important. This smooth-composition model can be 
applied also for calculating the electron transport proper- 
ties (thermal and electric conductivities) in the deep layers 
of the crust. We have calculated the static-lattice contri- 
bution for nonspherical nuclei at the very crust bottom, 
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at densities p <; 10 g cm -3 . We have analyzed (Sect. 
4) the neutrino emissivity as a function of density, tem- 
perature, and nuclear composition. We have obtained a 
simple analytic fit jic| ) to the bremsstrahlung emissivity 
for ground-state matter composed of spherical nuclei in 
the neutron-star crusts. We expect the fit to give reliable 
values of the neutrino bremsstrahlung emissivity for den- 
sities p < 10 14 g cm -3 , where the nuclei are certainly 
spherical, and to give a reliable upper limit of the emis- 
sivity for 10 14 g cm~ 3 < p < 1.4 x 10 14 g cm -3 , where the 
nuclei become nonspherical. 

Our results can be used, for instance, to study cooling 
of isolated neutron stars. A not too old star (of age lower 
than about 10 5 - 10 6 yr) cools mai nly vi a neutrino emis- 
sion from its interior (e.g., Pethick 1992 ). The main con- 
tribution to the neutrino luminosity comes usually from 
the stellar core. However, neutrino emission from the crust 
can also be important. In young neutron stars it plays a 
significant role in thermal relaxation of the stellar interi- 
ors. It can also be important in rather old neutron stars 
during the transition from the neutrino cooling stage to 
the photon one. The neutrino luminosity of the stellar 
core decays generally somewhat faster than the luminos- 
ity of the crust, and the crustal luminosity survives longer. 
Moreover, the neutrino emission from the crust can dom- 
inate the emission from the core in the neutron stars with 
highly supcrfluid cores and/or in stars with a stiff equation 
of state. In the latter case, the crust can be quite massive 
and its neutrino luminosity can be substantial. Finally, 
low-mass neutron stars always possess relatively massive 
crusts, and their neutrino luminosities can be mainly de- 
termined by their crusts. 

A number of problems remain to be solved in connec- 
tion with the calculation of neutrino pair bremsstrahlung 
in the crust. One of these is the influence on the phonon 
spectrum of neutrons outside nuclei when matter is made 
up of spherical nuclei. Intuitively one would expect the 
neutron liquid to make the effective mass of nuclei larger, 
thereby decreasing phonon frequencies. This in turn would 
increase the bremsstrahlung from phonon processes at 
least at low temperatures. Another is the nature of the 
collective excitations of the phases with non-spherical nu- 
clei. There are a number of these that have low frequencies, 
and these will give a "phonon contribution", while at the 
same time affecting the static lattice contribution. 
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